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We demonstrate a robust experimental method for determining the depth of individual shallow 
Nitrogen-Vacancy (NV) centers in diamond with ~ 1 nm uncertainty. We use a confocal micro¬ 
scope to observe single NV centers and detect the proton nuclear magnetic resonance (NMR) signal 
produced by objective immersion oil, which has well understood nuclear spin properties, on the di¬ 
amond surface. We determine the NV center depth by analyzing the NV NMR data using a model 
that describes the interaction of a single NV center with the statistically-polarized proton spin bath. 

We repeat this procedure for a large number of individual, shallow NV centers and compare the 
resulting NV depths to the mean value expected from simulations of the ion implantation process 
used to create the NV centers, with reasonable agreement. 


I. INTRODUCTION 

The Nitrogen-Vacancy (NV) center in diamond is a 
leading platform for wide-ranging applications in sensing, 
imaging, and quantum information processin^iH^. Key 
enabling properties of NV centers include excep tion ally 
long electronic spin coherence times (T 2 > 100/rs)pl and 
optical polarization and readout of the spin state (Fig. 

in an atomic sized defect within the diamond crystal 
under ambient conditions. 

Shallow NV centers within several nanometers of the 
diamond surface are especially useful for applications 
that rely on the strong dipolar coupling afforded by 
bringing the NV spin into close proximity to an exter¬ 
nal spin of interest. For example, quantum sensing 
and computinj^ schemes in which NV centers are em¬ 
ployed to control and read out the states of nuclear 
spins in samples tethered to the diamond surface re¬ 
quire minimal separation between the NV and nuclear 
spins for strong coupling. In magnetic sensing applica¬ 
tions, shallow NV centers with few nanometer separation 
from the magnetic field source have significant advan¬ 
tages over deeper NV centers and other magnetometers 
(e.g., SQUIDs) with much larger stand-off distances. Due 
to their close proximity to the sample, shallow NV centers 
(i) experience a larger magnetic held (i.e., dipolar helds 
fall off as l/r^) and (ii) enable spatial resolution on a 
length-scale c ompa rable to the stand-off distance, e.g., 
using scanninJ^ESI, super-resolution opticaP^^, or Fourier 
imagin^i^ techniques. In particular, shallow NV centers 
have recently been used for nuclear magnetic resonance 
(NMR) spectroscopy and magnetic resonance imaging of 


nanoscale sampleJ^^HIIl including single proton NMR and 

mrM 

Such applications of shallow NV centers depend cru¬ 
cially on accurate determination of the NV center depth, 
with uncertainty ~ 1 nm. Shallow NV centers are most 
commonly formed via nitrogen ion implantation, with the 
NV center depth estimated using the Stopping and Range 
of Ions in Matter (SRIM) Monte-Carlo simulatioil^. 
However, these estimates are statistical and thus do not 
provide the depth of any individual NV center. Further¬ 
more, the simulations do not take into account crystal¬ 
lographic effects such as ion channeling, leading to un¬ 
derestimation of the NV depth by as much as a factor 
of twcJi^. NV depth has also been estimated using sec¬ 
ondary ion mass spectroscopy (SIMS) of nitrogen ions 
after implantation. Unfortunately, SIMS has a minimum 
detection threshold (^ 3 x 10^"^ ^^N/cm^) and cannot be 
used to estimate individual NV center depthJi^. 

Recently, the depth of individual NV centers has been 
experimentally determined using two techniques requir¬ 
ing highly-specialized and delicate apparatus. The first 
technique takes advantage of Forster Resonance Energy 
Transfer (FRET), determining NV depth by observing 
the coupling of single NV centers and a sheet of graphene 
brought in close proximity with the diamond surface. 
Measuring the NV fluorescence intensity as a function 
of separation between the graphene and diamond surface 
until the two are in contact and fitting the data with 
a theoretical model, NV depth can be determined with 
sub-nanometer uncertainty^. In the second technique, 
a single shallow NV is employed to image, with ^ 1 nm 
vertical resolution, dark electron spins assumed to be lo- 
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(^) (b) Nuclear spin 



FIG. 1: NV NMR Experiment, (a) NV electronic en¬ 
ergy level structure, (b) A confocal microscope addresses a 
single shallow NV center, which detects NMR signals from 
a few-nanometer region of sample on the diamond surface. 
Due to dipolar coupling, a shallow NV center (left) experi¬ 
ences a signihcantly stronger magnetic field from a smaller 
nuclear spin sample volume than a deep NV center experi¬ 
ences (right). The strength of the magnetic field at the NV 
center is indicated by the opacity of the nuclear spin sample, 
and the dashed lines qualitatively illustrate the volume of nu¬ 
clear spin sample that contributes most of the NMR signal, 
(c) Larmor precessing statistically-polarized nuclear spins in 
the sample produce an effective AC magnetic field (green) 
that is detected by the NV sensor in a frequency-selective 
manner using an XYSfc pulse sequence. 

cated at the surface of the diamond sample. The dark 
spin imaging resolution and consequently the uncertainty 
in NV depth determination is ultimately limited by the 
applied magnetic field gradient, the mechanical stability 
of the apparatus, and the T 2 of the dark spirP^. 

In this paper, we present a robust method for ex¬ 
tracting individual NV center depth with ~ I nm un¬ 
certainty that can be easily performed with a scanning 
confocal microscope. We derive and analyze a model that 
describes the interaction of a single shallow NV center 
with a statistically-polarized nuclear spin bath, such as 
a proton-containing sample on the diamond surface, and 
discuss the conditions of validity of this model. Fitting 
the single-NV-measured proton NMR signal produced by 
microscope objective immersion oil, which has well un¬ 
derstood nuclear spin properties, to the model expres¬ 
sion, we determine depths for a large number of individ¬ 
ual shallow NV centers and compare the measured depths 
with those expected from SRIM simulations. Finally, we 
discuss further application of this model to perform char¬ 
acterization of both NV centers as well as unknown nu¬ 
clear spin samples on the diamond surface. Note that 
the experiments, model, and analysis presented here are 
a more detailed treatment of similar approaches to de¬ 
termining NV depth outlined in Refs. IT1^I5I2()1 

II. METHODS 

In our experiments we study negatively-charged NV 
centers formed via low-energy, low-dosage nitrogen ion 
implantation and subsequent annealing (see details in 


Sec. Ill and Table |I]), such that individual NV centers 
can be interrogated with a confocal microscope. To de¬ 
termine the depth of an individual NV center, we apply 
immersion oil to the diamond surface and measure the 
variance of the fluctuating NMR magnetic field at the 
NV center using a dynamical decoupling pulse sequence. 
The NMR magnetic field is created by a statistically- 
polarized subset of the proximal protons in the immersion 
oil, as shown in Fig. [^. The protons undergo Larmor 
precession with a frequency determined by the applied 
static magnetic field (I50-I600 G), but with a phase and 
amplitude that varies with every repetition of the pulse 
sequence. Although the net magnetization of the proton 
spin ensemble over the timescale of the entire experiment 
is negligible at the temperature and static fields applied 
in this work, the variance is nonzero and is proportional 
to the density of the proton bath. 

We use an XYSfc pulse sequence, shown in Fig. [J, 
to measure individual Fourier components of the NMR 
magnetic signal. We first optically pump the NV center 
electronic spin into the TOs = 0 magnetic sublevel and cre¬ 
ate a coherent superposition of the m-s = 0 and m-s = I 
sublevels using a microwave (MW) 7r/2-pulse. The NV 
spin then undergoes periodic intervals of free evolution 
and 180° phase flips driven by resonant MW pulses, af¬ 
ter which a final MW 7r/2-pulse converts the accumulated 
phase into an NV spin state population difference. The 
NV spin free evolution is governed by the time-dependent 
component of the total external magnetic field, which in¬ 
cludes contributions from the proton NMR signal pro¬ 
duced by the immersion oil on the diamond. The net 
accumulated NV spin phase is only appreciable when the 
evolution time r is close to half the proton Larmor pe¬ 
riod. 

The accumulated NV spin phase is measured by two 
consecutive near-identical experiments that project the 
final NV spin state first onto the TOs = 0 state (re¬ 
sulting in a measurement of NV fluorescence Fq) and 
then onto the rris = 1 state (resulting in a measure¬ 
ment of NV fluorescence Fi), with appropriate choice of 
the final 7r/2-pulse phase. In order to remove common¬ 
mode noise from laser fluctuations, the two fluorescence 
signals are normalized to give the signal contrast S = 

[(Fo-Fi)/(Fo + i^i)]. 

Measuring the signal contrast over a range of free evo¬ 
lution times T results in slowly decreasing signal contrast 
for larger r, due to NV spin decoherence, and a narrower 
dip in contrast for specific values of r, caused by the nu¬ 
clear spin Larmor precession. The background decoher¬ 
ence can be fit to an exponential function and normalized 
out, leaving the normalized contrast C(t) with only the 
narrower NMR-induced dip (shown in detail in the ap¬ 
pendix). The shape of this dip, described by Equation]^ 
is determined by the magnetic field fluctuations produced 
by the dense ensemble of nuclear spins in the immersion 
oil on the diamond surface, as well as by the filter func¬ 
tion corresponding to the XYSfc dynamical decoupling 















3 


pulse sequence: 


C{t) « exp 





( 1 ) 


(An in-depth derivation is presented in the appendix.) 
Here 7 e « 1.76 x 10^^ rad/s/T is the electron gyro- 
magnetic ratio, B^ms is the RMS magnetic field signal 
produced at the Larmor frequency by the nuclear spins, 
K{Nt) is a functional which depends on the pulse se¬ 
quence and the nuclear spin coherence time, and N is 
the number of 7r-pulses, which are separated by the NV 
spin free precession time r. As shown in the appendix, 
for the simplest case of a semi-infinite layer of a homoge¬ 
neous nuclear-spin-containing sample on the most com¬ 
monly used {100}-oriented diamond surface, Hrms is re¬ 
lated to the NV depth d^v below the diamond surface 
by 
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where p is the nuclear spin number density and 7 „ is 
the nuclear spin gyromagnetic ratio (for protons 7 „ « 
2.68 X 10® rad/s/T). More general cases of arbitrary 
nuclear spin quantum number and other diamond sur¬ 
face orientations can be calculated as described in the 
appendix. If the nuclear spin dephasing time (T^n) is as¬ 
sumed to be infinite, then the functional K{Nt) is given 

by 


K{Nt) 


(Vr)^sinc^ 


Nt ( 



( 3 ) 


where w/, is the nuclear Larmor frequency. However, 
spectral broadening of the NMR signal due to diffusion 
or a finite dephasing time can also be included as shown 
in the appendix, in which case, the functional K{Nt) is 
given by 



For a sample with well-known nuclear spin number den¬ 
sity p (e.g., p = 68 ± 5 nm“® for the Nikon Type NF 
immersion oil employed in this work, measured using a 
Varian Unity InovaSOOC NMR system), the only free pa¬ 
rameters in the fit expression are the NV depth the 
Larmor frequency and the nuclear spin dephasing 
time The confidence with which each of these pa¬ 
rameters can be extracted from a fit of Equation to NV 
NMR data is strongly dependent on both the probed NV 
center properties and the applied pulse sequence. 

In the limit of infinite T|„, the strength of the NMR 
signal dip is entirely determined by the NV depth and 
the measurement pulse sequence duration T = Nt, vary¬ 
ing inversely with the former and directly with the latter. 
That is, for a fixed pulse sequence duration, shallower NV 
centers produce stronger NMR signal dips while deeper 
NV centers produce weaker NMR signal dips. As a re¬ 
sult, pulse sequences with longer durations are necessary 
to acquire a strong enough NMR signal dip to confidently 
extract a depth estimate from a deeper NV center. On 
a related note, the infinite limit is only valid when 
the pulse sequence duration is significantly shorter than 
T|„; for sufficiently long pulse sequence duration, the NV 
detection bandwidth becomes narrow enough that the 
broadening of the NMR signal dip due to nuclear diffu¬ 


sion and spin dephasing can be observed and can be 
extracted using the form of the functional K{Nt) given 
by Eq. The pulse sequence duration is eventually 
limited by the coherence time T 2 of the NV spin, however, 
which places upper bounds on the depth of NV centers 
and T 2 „ of nuclear spin samples that can be extracted 
with this analysis. Recent work indicates a strong de¬ 
pendence of the NV T 2 coherence time on the NV depth 
for shallow NV centers.^ Assuming a typical value of 
T 2 ~ 1 ms found in deep NV centers and standard opti¬ 
cal collection efficiencies (< 10%) we estimate that NV 
depths up to 300 nm below the diamond surface can be 
measured using the present method. 


III. RESULTS 

We performed measurements on 36 NV centers across 
3 diamond samples, each synthesized via chemical vapor 
deposition (Element Six) and isotopically engineered to 
contain 99.999% Sample A was implanted with 3- 
keV ^®N+ ions at a dose of 1 x 10® cm“®; Sample B was 
implanted with 2-keV ^^N®" ions at a dose of 1 x 10® 
cm“^; and Sample C was implanted with 2.5-keV ^^N+ 
ions with measurements taken in a region of 2D NV den- 
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FIG. 2: Example NV NMR proton spectra. For all spectra, diamond sample and NV pulse sequence, and applied 
static magnetic field are given in the bold inset label or in the symbol key, and the extracted NV depths are given in the symbol 
key. (a) NV NMR proton spectra data (black dots) measured with an XY064 pulse sequence at 197 G static field, analyzed 
assuming finite T 2 „ (red solid curve) and infinite r|„ (blue dashed curve). Both analyses fit the data well, with consistent NV 
depth values, (b) Proton NMR spectra measured with another NV center using different pulse sequences. The NV depths 
extracted from finite T^n fits (solid curves) are in reasonable agreement for all measurements, (c) Proton NMR spectra and 
finite T2n fits (solid curves) for two NV centers determined to have different depths under the same experimental conditions. 
The observed signal contrast dips vary strongly with NV depth, (d-f) Proton NMR spectra measured with the same NV center 
at different static field strengths and using different pulse sequences. Finite fits (solid curves) yield consistent NV depths 
for all experimental conditions. 


sity ^ 8 X 10^ cm“^. We employed a custom-built scan¬ 
ning confocal microscope to address single NV centers 
in each sample and fit the measured proton NMR signal 
from immersion oil on the diamond surface to Equation 
[2 in order to extract depth values for each NV center. 
A compilation of the measured properties of all the NV 
centers and diamond samples is given in Table |l] Pro¬ 
ton spins in immersion oil have an expected ~ ^0 /rs 
(corresponding to a linewidth ~ 5 kHz, see appendix for 
details) which is a longer nuclear than can be ex¬ 
tracted with the shallow NV centers used in the present 
work. Indeed, analysis of the measured NMR spectra 
data assuming infinite (Eq. and finite (Eq. 
generally give good agreement both in fits to the data 
and in NV depth extracted (Eig. However, since 

the infinite condition does not hold strictly true for 
every measurement, we performed all analyses using the 
general case of finite nuclear except where explicitly 
noted. 

Eigure shows typical measured proton NMR data 
from several representative NV centers in Sample A. The 
solid curves correspond to the best-fits of the model func¬ 
tion to the data, from which NV depth estimates are 
extracted. We find that the contrast dip positions are 


in good agreement with those expected for the magnetic 
fields measured from the NV resonance frequencies, i.e., 
dips occur at r = tt/lul. Furthermore, we find that the 
fit expression yields consistent NV depth values even un¬ 
der different experimental conditions. For example, in 
Figure [^, several measurements with different numbers 
of pulses were performed on the same NV center at the 
same static magnetic field. Fitting to each NMR spec¬ 
trum independently, we extracted NV depth values that 
were in reasonable agreement with each other. Figures 
id -f) show measurements and analyses of another NV 
center for which both the number of pulses and the static 
magnetic field were varied. Again, for all experimental 
conditions, the NV depth values extracted from the mea¬ 
surements are comparable to within their error bars. Fig¬ 
ure [^c) shows proton NMR data from two different NV 
centers measured with the same pulse sequence under the 
same experimental conditions (within the same diamond 
sample at the same static magnetic field) to illustrate 
the profound effect an NV center’s depth can have on 
its sensitivity to NMR signals from nuclear spins at the 
diamond surface. 

Finally, we compared the distribution of NV depth val¬ 
ues extracted from diamonds with different nitrogen im- 
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FIG. 3: Histogram of measured NV depths in two di¬ 
amond samples. Estimated depths of (a) 11 NV centers in 
diamond Sample A, implanted with 3.0-keV ions, and (b) 
13 NV centers in diamond Sample C, implanted with 2.5-keV 
^^N ions. 


plantation energies. Figure shows histograms of the 
estimated depths for 11 NV centers in Sample A, which 
had been implanted with 3.0-keV ions and 13 NV 
centers in Sample C, which had been implanted with 2.5- 
keV ^^N ions (see also Table [^. We found that the 3.0- 
keV implanted NV centers had a mean depth of 10.5 nm, 
with 2.8 nm standard deviation, and that the 2.5-keV im¬ 
planted NV centers had a slightly shallower mean depth 
of 8.5 nm, with 2.8 nm standard deviation. In contrast, 
SRIM simulations predict a mean depth of (5.2 ± 2.1) 
nm for 3.0-keV ion implantation and a mean depth 
of (4.5 ± 1.9) nm for 2.5-keV; thus our measurements 
of NV depth are consistent with previous estimates that 
SRIM underestimates NV depth by as much as a factor 
of two.^ However, it is important to note that the SRIM 
software estimates the distribution of implanted nitrogen 
ions whereas the NV NMR analysis estimates the depths 
of NV centers, which may have depth-dependent factors 
limiting their formation in diamond beyond the distribu¬ 
tion of implanted nitrogen impurities. Furthermore, in 
addition to the NV centers whose extracted depths are 
represented in Figure in all diamond samples we ob¬ 
served that a fraction of the optically observed NV cen¬ 
ters (e.g., roughly 1/2 in Sample C) had optical and/or 
spin properties that were too unstable for any detailed 
measurements to be performed on them. These unstable 
optical and/or spin properties are likely symptomatic of 
very shallow NV centers whose depths cannot therefore 
be measured with the NMR technique presented in this 
paper. While this behavior may indicate a bias in the NV 
depth statistics extracted using this analysis technique, 
it also illustrates how this analysis may be applied to¬ 
wards determining how close to the diamond surface NV 
centers’ optical and spin properties remain stable enough 
for sensitive spin measurements and furthermore provides 
an avenue for studying how surface treatments and pro¬ 
cessing can be used to stabilize very shallow NV centers. 
Both are topics of great importance in sensing, imaging, 
and quantum information applications that rely on shal¬ 
low NV centers. 


TABLE I: Summary of the depths determined from 36 NV 
centers in 3 diamond samples under a range of external static 
field magnitudes Bo and number of vr-pulses N used in the 
XY8A: measurement protocol. Sample A was implanted with 
3.0-keV ^®N ions; Sample B was implanted with 2.0-keV ^®N 
ions; and Sample C was implanted with 2.5-keV ^'’N ions. In 
Samples A and C, measurements were performed on a ran¬ 
dom collection of NV centers such that the determined depth 
values reflect the NV depth distribution. In Sample B, mea¬ 
surements at 1609 G were performed only on NV centers that 
showed strong proton NMR signals for short averaging times; 
consequently these measurements are weighted towards shal¬ 
lower NV centers and do not accurately reflect the NV depth 
distribution. 


Sample 

NV# 

Bo (G) 

TT-pulses 

NV depth (nm) 

A 

001 

197 

32 

10.4(7) 

A 

002 

197 

64 

13.2(3) 

A 

005 

197 

64 

14.8(3) 

A 

006 

197 

16, 32, 64 

8.5(4) 

A 

007 

197 

32, 64 

9.0(4) 

A 

008 

197 

64, 256 

15.3(3) 

A 

010 

197, 1580 

16, 32, 64, 508 

8.9(5) 

A 

012 

197 

32 

8.3(3) 

A 

104 

150 

16 

6.4(2) 

A 

110 

150 

64 

10.7(4) 

A 

111 

150 

64 

10.0(2) 

B 

009 

206 

64 

10.7(7) 

B 

022 

159 

32, 64, 96, 128 

9.7(6) 

B 

100 

206 

32 

11(2) 

B 

112 

1609 

60 

6.2(6) 

B 

115 

1609 

124 

7.7(3) 

B 

116 

1609 

124 

5.2(2) 

B 

118 

1609 

124 

6.5(3) 

B 

119 

1609 

124 

4.8(2) 

B 

120 

1609 

124 

4.8(2) 

B 

121 

1609 

124 

5.6(3) 

B 

122 

1609 

124 

5.0(2) 

B 

123 

1609 

124 

7.3(3) 

G 

009 

156 

16, 32, 64, 96 

8(1) 

G 

014 

156 

64 

13.3(9) 

G 

025 

156 

64 

9.4(5) 

G 

030 

156 

16 

4.9(4) 

G 

056 

156 

8, 16 

4.7(2) 

G 

075 

156 

64 

7.4(2) 

G 

090 

156 

64, 96, 128 

7.5(5) 

G 

093 

156 

64, 128 

9.4(6) 

G 

098 

156 

64, 96, 128 

12(1) 

G 

107 

156 

64 

8.6(4) 

G 

111 

156 

16, 32 

4.6(6) 

G 

116 

156 

64 

9.7(6) 

G 

125 

156 

64, 128 

11(1) 


IV. DISCUSSION 

Our robust NMR technique for determining the depth 
of shallow NV centers also enables detailed investigations 
of the effect of NV depth on other NV center properties. 
In particular, NV spin properties such as dephasing time 
T 2 , coherence time T 2 , and relaxation time Ti may be 
characterized as a function of depth; furthermore, NV 
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spectroscopic techniques may be applied to probe the 
local spin environment close to the diamond surfac^^. 
Since magnetic sensing and quantum information appli¬ 
cations that employ shallow NV centers also require long 
NV spin coherence times, better understanding and con¬ 
trol of NV spin properties and the spin environment as a 
function of NV depth are key challenges. 

In the present work, we applied our technique to deter¬ 
mine NV center depth using a well-known nuclear sample. 
However, once an NV center’s depth is determined, this 
information can be combined with our model to perform 
NV NMR studies of unknown nuclear samples. Also, as 
discussed in Section III, applying appropriate pulse se¬ 
quences allows for the extraction of the nuclear spin 
which can be used to study nuclear spin interactions and 
diffusion in the sample. Furthermore, by probing an un¬ 
known nuclear sample using multiple NV centers of dif¬ 
fering depths, information about the nuclear spin distri¬ 
bution as a function of sample depth may be extracted.!^ 
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FIG. A.l: (a) Example NV signal contrast Sir) data (circles) 
measured by applying an XY064 pulse sequence on NVA006 
(Sample A). The decay due to NV spin decoherence is fit to 
a stretched exponential function (line), excluding the data 
which makes up the narrow NMR dip (open circles), (b) 
Normalized contrast C(r) data isolates the NV NMR signal. 


Appendix B: NV NMR Lineshape 
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Appendix A: NV Spin Decoherence Normalization 


As described in the main text, two NV“ spin-state- 
dependent fluorescence measurements Fq{t) and Fi{t) 
are acquired from consecutive, near-identical but inde¬ 
pendent dynamical decoupling experiments, each with 
TT-pulses spaced by time r. For Fb(r), the final 7r/2- 
pulse projects the NV spin coherence onto the |0) state, 
whereas for F’i(r) the pulse phase is reversed to project 
the coherence onto | ± 1). This procedure removes 
common-mode noise from laser fluctuations occurring on 
timescales > r. The fluorescence signals are described as 
a signal contrast, S'(r), of the form: 


S[t) = 


Fo{t) - Fi{t) 
Fair) +Fi{t)' 


(Al) 


The signal contrast effectively measures the projection of 
the NV spin coherence after the pulse sequence onto the 
coherence at the beginning of the sequence. Measuring 
S over a range of free evolution times r yields a slow de¬ 
cay due to NV spin decoherence and a narrow dip due 
to nuclear spin Larmor precession. The background NV 
spin decoherence can be fit to a stretched exponential 
function, excluding the data points which make up the 
narrow dip corresponding to the NMR signal, as shown 
in Fig. A.l[ a). Dividing by this exponential fit func¬ 
tion yields a normalized contrast C'(r) which isolates the 
NMR signal in the NV measurement, as shown in Fig. 
^b). 


In this appendix, we present a derivation for the sig¬ 
nal expected from an NV NMR measurement made with 
a dynamical decoupling sequence. We adopt the non¬ 
unitary Fourier transform in angular frequency units, 
such that the Fourier transform pair for f(t) is defined 

aPl 

fit)=^-\fiuj)) = - fiuj)F^^dw, 

/ + 00 

/(f)e-“*df. (Bl) 

-OO 

With the previous expression, Parseval’s theorem reads 
as: 



f{t)g*{t)dt 
\fit)\^dt = 


1 

27r 



f{u})g*{u})duj 


1 

27r 



l/(w)pdw. 


(B2) 


and the expressions for the Dirac delta and convolution 
functions are: 


Situ -uj') = — / 

27’' J-oo 

*g) = /(w)g(a;). 


(B3) 


1. Signal from a Dynamical Decoupling Sequence 


During the dynamical decoupling measurement se¬ 
quence, the NV spin coherence accumulates some phase 
A0 (t) due to evolution in the presence of magnetic fields. 
In this work, the magnetic field of interest is the NMR 
signal from statistically-polarized spins in the sample on 
the diamond surface. After normalizing out contributions 
due to background NV spin decoherence (see Appendix 
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A), the contrast is related to the accumulated phase by: 


C{t) = (cos(A(/)(r))). (B4) 


The brackets around cos(A(/)(r)) indicate that a typical 
fluorescence measurement is an average over many re¬ 
peated, nominally identical dynamical decoupling exper¬ 
iments. If the accumulated phase A(/)(r) follows a nor¬ 
mal distribution centered at zero with variance (A(/)^(t)) 
as will typically be the case for an NMR signal from a 
statistically-polarized nanoscale sample, then the aver¬ 
age over the cosine can be converted to an exponential 
function of the variance using the relationshijPSI 



FIG. B.l: The dynamical decoupling sequence, induced by 
resonant MW pulses with phases as labeled, defines a function 
g{t) describing the direction of NV spin precession in response 
to a magnetic signal Bz{t). 


/ C30 

f{x)p{x)dx, (B5) 

-OO 


where p{x) is the probability distribution function for 
random variable X. Applying the integral of Eq. (B5) to 
Eq. (|B4[) yields: 


Then we can write: 


C'(r) = exp(-(A(/)2(r))/2). 


(B6) 


Phase is accumulated during the dynamical decoupling 
sequence as the NV electronic spins Larmor process in the 
presence of a magnetic field signal Bz{t), where z is the 
NV quantization axis. (The NV spin Larmor precession 
from the static background field Bq is removed by work¬ 
ing in the rotating reference frame). The sign of phase 
accumulation (i.e., positive or negative phase accumula¬ 
tion) is reversed by each 7r-pulse of the sequence, and can 
be represented over time as a function g{t), as shown in 
Fig. B.l The total phase accumulated at the end of the 


sequence is then: 


/ -Too ^-1-00 
-c 


' —OO —OO 

/ +00 /*+oo 

/ 

-OO —OO 

/ +00 />+oo 

/ 

-OO j —OO 

/ + 00 

JUt')9it')dt' 

-OO 

/*+oo 

-^i,2(w)5(w)dw 


SB{T)g{t')g{T -h t')dTdt' 

SB{T)g{T -h t')dTg{t')dt' 
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27r 


(BIO) 


where in the last line of the previous expression we have 
used Parseval’s theorem. Since the term Jf 2(^0 noth¬ 
ing but a convolution, one can easily conclude that: 


^ ^2 /--l-oo 

{^4> ij)) = ^ J SBixi)g{-uj)g{uj)du 


7^ 

27T 


— OO 
2 /•+00 


SBiuj)\g{ui)fduj. 


(Bll) 


A(^(r) = 7e 


^+oo 


g{t)Bzit)dt, 


(B7) 


where 7 e is the gyromagnetic ratio for the NV electronic 
spin (in units of rad/s). The accumulated phase vari¬ 
ance can be expressed in terms of a correlation function 
between measurements across times t and t': 

/ -l-oo l•+oo 

g{t)Bz{t)dt / g{t')Bz{t')dt'). 

-OO j —OO 

(B8) 


The quantity Sb{pj) represents the spectral density of the 
effective NV spin phase noise resulting from the magnetic 
field Bz (t) and manipulation of the NV spin by repeated 
dynamical decoupling sequences; it can be computed as 
follows: 


Sb{uj) = {\Bz{u;)\^) 

r+oo 

= / {BziO)Bz{t'))e-^‘^^'dt'. 


(B12) 


2. Application to NMR Signals 


We now assume temporal translational invariance for the 
local and time-dependent field correlator: 


Correlation Functions 


We consider the NMR magnetic signal Bz{t) originat- 
(B9) ing from nuclear spins on the surface of the diamond 


{Bz{t)Bz{t'))=SB{t-t'). 






















and in the vicinity of a shallow NV center (see Fig. 


B.2). The statistically-polarized nuclear spin ensemble 


produces fluctuations in (t). For an ensemble of point 
dipoles, Bz{t) at the NV center can be written as: 


Bz{t) = ^ D, [iuWzIiit) + 

3 

+ (3w>i-l)/i (t)] , (B13) 

where the NV is coupled to many nuclear spins j at posi¬ 


tions given by a distance Vj and a unit vector (which 
can be written in terms of of its coordinates ui.,u^y,ui). 
The coupling factor is Dj = (/io/i 7 „)/( 47 rr|), where 7 „ 
is the gyromagnetic ratio of the nuclei and rj is the dis¬ 
tance between the NV center and nuclear spin j. Terms 
^x,y,z represent the operator projection of nuclear spin j 
along the x, y, and z axes. 


Using Eq. (B13), the time-dependent correlator for the 
NMR magnetic field can be expressed as: 


{Bz{0)Bz{t)) = [3u>i/^(0) + 3«/^(0) + (3wi< - l)Ji(O)] 

3 

Y^D^Vi) [3«/l(t) + 3u>*7;(t) + (3u>i - l)Ilit)]). 


(B14) 


For an ensemble of nuclear spins that do not interact 
with each other, time-dependent correlators can be de¬ 
fined for every spin’s operator projection along each of 
its axes: 

iimipit)) = (B15) 

The function represents the local nuclear spin-spin 
correlation function. By treating the nuclear spins as 
paramagnetic, the correlations between different nuclear 
sites are identically zero. Note that the correlator is a 
function of the nuclear spin’s total spin quantum num¬ 
ber I as well as the temperature T and the applied field 
Bq (which determines the Larmor frequency of the nu¬ 
clei). In the simple case in which the external magnetic 
field for the nuclei is applied along the NV axis one can 
write ^ i.e., behavior in the transverse plane 

is independent of the relative phase between the nuclear 
spin and the NV. Moreover, all nuclear spins of the same 
species have the same correlator, and so the index j is 
dropped for Then 

{Bzim^it)) [ 9 /"’" {{nWzf + {u^zf) 

3 

{suiui-if. 

(B16) 


Then in the spectral representation 


r+oo 


/“’“(/,T,a;) = ^(/“’“(t)) = 


(/„(t)J„( 0 ))e-“‘dt 


27T 




where Z is the spin partition function and „ are the 
energies of nuclear spins m,n!^ In the high temperature 
limit where En <C ksT, the eigenstates are equally pop¬ 
ulated, and 




27r 


Y fz\Ia\mz) 

n,m 


(B19) 


We now make use of the definitions for the z and x spin 
projections: 


Iz = {nz\Iz\mz) = mz{nz\mz) 

Ix = {nz\ ^ \mz), (B20) 


where 

I^\I,mz) = y/pl + 1 ) - mzimz ± ± 1 ). 

(B21) 


Assuming that the energy of the nuclear spin state \mz) 
is HuJm^rUz, the transverse f^’^, and longitudinal 
spin-spin correlation functions have their natural 
expression in frequency-space with the definition in Eq. 
(Bl|. The relevant spin projections la for each nucleus 


are found using their respective operators: 


la = 




(B17) 


Then the longitudinal correlator is 

|(m.|J.|m.)|"<5(.;). (B22) 


The correlator ( B22[ ) can be computed by noting that 
a Curie-Weiss prefactor appears due to the relation 
Ylz ™ 0 /Tr(l) = /(/-l-l)/3. Because the longitudinal cor¬ 
relator is centered at zero energy, it will not contribute to 
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FIG. B.2: An NV center at depth d below the diamond surface 
on which resides a sample containing an ensemble of nuclear 
spins, each with spin vector Ij and position ul.,Uy,ui. The 
NV axis, and the axis for magnetic quantization, is at angle 
a with respect to the vector normal to the diamond surface. 
For purposes of integration across the sample, the spherical 
coordinates r,9,cj> are used. The external magnetic field Bq 
is assumed to be aligned with the N-V symmetry axis. 


the final integral (Bll I as long as g{u! = 0, r, iV) = 0 (i.e., 
the dynamical decoupling pulse sequence is not sensitive 
to DC fields). The transverse correlator is 


27r 






- a; 


(B23) 

which is non-zero only when mz^n^ are adjacent en¬ 
ergy levels. For the case of spin-1/2 nuclei (I = 1/2), 
where the nuclear spins precess at Larmor frequency 
wl = 7n.So, we evaluate ( B23| ) as: 

O-TT 

= 1/2, a;) = — - ujl) + <5(a; + w^)). (B24) 


The two contributions in (B24) represent the Stokes and 
anti-Stokes lines, equal in the limit T oo.l^ 

The expression for magnetic field correlation is now 

{Bzm.{t)) = + [uiuif] , 

3 

(B25) 

with given by Eq. B24 By writing 1 — (u^)^ = 
+ (w^)^j the geometry-dependent terms can be col¬ 
lected into one factor: 




(B26) 


which we evaluate in the following section. 


b. Calculation of the Geometrical Factor 

For liquid samples such as immersion oil in which nu¬ 
clear locations vary on a time scale short compared with 
the dynamical decoupling sequence length, one can as¬ 
sume a sample of nuclear density p continuously dis¬ 
tributed on the diamond surface. Then the summation 


of the geometrical factor (B26) can be converted to the 
integral: 


V = p dV 


\ 47r 



(B27) 


We evaluated the integral F using spherical coordinates 
with the conventions of Fig. |B.2| The polar angle ori¬ 
gin 0 = 0 is defined to be orthogonal to the surface of 
the diamond, while (f> is the azimuthal angle with arbi¬ 
trary origin. The NV axis z points along a direction 
z = [sin(a) cos(/3) , sin(a ) sin(/3), cos(q)]. The projection 
Uz needed for Eq. ( B27| ) will in general depend on all four 
angles just introduced. In particular, Uz = z • u^, where 
Ur = [sin(0) cos((/), sin(0) sin((/), cos(0)]. 

The integral for F is then 


n27T /*7r/2 pOO 


F = 


iuz)^(i-{uzr) 


to Jo 


' d]S!V / cos(0) 


sin(0)drd0d(/), 

(B28) 


where (Inv is the NV depth below the diamond surface. 
The sample height is assumed to be semi-infinite, thereby 
allowing integration of the radial component from the 
diamond surface to infinity. Other sample geometries 
can be accomodated with the proper integral limits and 
choice of coordinate system (i.e., spherical, cylindrical, 
etc.). Evaluating the integral produces a simple expres¬ 
sion for r(fiAry): 


r(dw) = P 


^ / TT [8 — 3sin^(a)] \ 

y 288d^Y j 


(B29) 


The expression is maximal when a = 0, where F(djvv) = 
TT/{36d%Y) However, in most diamond samples, the nor¬ 
mal to the surface is aligned along the [100] crystal di¬ 
rection, so that a = 54.7°. At this angle, r(djvy) = 
57r/(216ci^y). With the correlation functions and geo¬ 
metric factors now evaluated, the spectral density can be 
written as: 


^B(a;) = (|S.(/ = 1/2 ,cc)|2) 

Qtt 

= T{dNv)^ - ujl) + S{uj + ujl)) ■ (B30) 

The spectral density can be related to the magnetic field 
variance from the NMR signal by: 

<S'b(w) = {d{uj — LOif) -\- -\- ojl)) , (B31) 


where 


B, 


RMS — 7r(dArv) 


_ / Poh'^n y / ^ [8 - 3 sin'^(a)] 


V 47 r 


128d^y 


(B32) 
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For NV centers oriented at a = 54.7° this simplifies to: 


1d2 

-^RMS 


f /i0^7n \ ' 


f 57r 


\96d^ 


NV 


(B33) 


If the nuclear spin sample on the diamond surface is semi¬ 
infinite laterally but not vertically, such as a thin layer 
between coordinates zi and Z 2 above the diamond sur¬ 


face, then Eq. B32 can be rewritten as: 


^RMS — P 


(PofilnV 3sin'‘(Q;)] 


V 47r 


1 


128 

1 


{dNV + {dNV + ^2)^, 


(B34) 


c. The Filter Function \g{uj,T)\‘^ 

To complete evaluat ion o f the accumulated NV spin 
phase variance integral (Bll) and thus the signal contrast 
Eq. ( |B6[ ), the filter function \g{uj, r)p must be determined 
for the dynamical decoupling sequence. For a CPMG or 


XY8 sequence with N 7r-pulses, such as that in Fig. Ic, 
we compute the Fourier transform: 


/ 2 ^ .Wx/ (2fc + l),r \ 


k— — c 


Sine 


Nt / {2k + l)7r 

—— UJ — - 

2 V r 


(B35) 


For most purposes, only the first-order terms in Eq. (B35) 


need to be retained. Additional terms contribute only to 
higher harmonics, which are not measured in this work. 
The expansion must include k = 0,-1 to be symmet¬ 
ric around iw. However, the integral over positive and 
negative frequencies will be equivalent to twice the in¬ 
tegral over positive frequencies as long as /cbP S> Hujl. 
If the nuclear spin dephasing time is assumed to be infi¬ 
nite, such that the nuclear spin signal can be described 
by delta functions, we can now obtain a final formula for 
the signal contrast in the 1=1/2 case, keeping terms 
k = 0,-1: 


C'(t) « exp 


~ ^RMS 


sme 


Nt 

~Y 


(-- 7 ) 


• Sine 




2 sine 


Nt 

( 


Nt / 7r\ 


~Y 


sine 







( 


(B36) 


The off-resonant terms contribute very weakly to the line- 
shape and can be ignored, resulting in an approximate 
formula: 

C{t) « exp -^'^lB\^^{NTfsm(? (ujl “ 

(B37) 


3. Nuclear spin dephasing time 

In the previous section, we assumed that the nuclear 
spin signal could be represented by a delta function, 
meaning that it has a dephasing time (T^n) much longer 
than the length of the NV dynamical decoupling se¬ 
quence. However, the effective nuclear spin linewidth is 
broadened due to both dephasing from spin-spin interac¬ 
tions and diffusion through the nanoscale NV interaction 
volume. In order to take these effects into account, we 


substitute the delta functions of Eq. (B241 with normal¬ 
ized Lorentzian functions such that: 


r^{i = 1 / 2 , 00 ) = 


2-k fl 


Tr 


2n 


8 \t{uo-ool)^ + {t;-^)^ 

1 rpif-l 

i- -t <2j^ 

+ V{~uo + uo^)^ + {t;-^^Y 


(B38) 


As before, we need to compute: 


= j ^’""{1,^) |5(W,T, V)|^dw^ . 

(B39) 

Once again, symmetry allows us to simplify the expres¬ 
sion using only the positive-frequency component if we 
multiply the expression by two, leading to: 
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C'(t) = exp 


— 1 


~;^7e'®RMS 




(iVr)^sinc 


Nt 

~Y 




(B40) 


It is evident that the integral is a convolution between a 
Lorentzian l(io) and a function ii){io) ^ sinc^(M). Using 
the convolution theorem, the integral can be solved by 
multiplying the respective Fourier transforms and then 
taking the inverse Fourier transform of the result. The 
Lorentzian component is 


1{U}) 


1 

TT (w 


Tl 


.-1 




(B41) 


Its Fourier transform is 


m = (, 


_ /g-tTin -itbJL 


(B42) 

where H{t) is the Heaviside step function. The sinc^(u) 


component is 


tpiuj) = (A^r)^sinc^ 


Nt 


(w) 


(B43) 


Notice that the frequency offset tt/t has been removed 
to simplify the Fourier transform. The Fourier transform 
is 


4'(t) = TT [(t — Nt) sgn(t — Nt) 

—2t sgn(t) + (t + Nt) sgn(t + iVr)] . (B44) 

Taking the inverse Fourier transform K(uj) = 
{L{t)'^{t)), and using the identity oj = tt/t for the 
filter function resonance condition, gives the expression: 


K{t) 





- 2r*„ (wi - ^ 



COS 

Nt (ojl - -) 

[ V T / 


L V r / J 


sin 


Nt {ujl 




1 + 




(B45) 


The final expression for signal contrast, including nuclear 
spin dephasing and again ignoring off-resonant terms in 
the filter function, is 


C{t) « exp ■ (B46) 

In practice, experimental determination of whether the 
nuclear spin is long or short relative to the length 
of the NV dynamical decoupling sequence can be carried 
out by checking the scaling of the observed contrast dip 
amplitude and width as a function of N and r. 


4. Pseudospin Derivation 

An alternative derivation of the signal contrast C'(t) 
can be obtained using the pseudospin formalisrrPI. The 
contrast is a product of the pseudo-spin signal from 
each nuclear spin j in the sample on the diamond surface: 

C'(r) = n5T (B47) 

3 


For a CPMG sequence (or XY8) with N pulses, the 
pseudo-spin signal for nuclear spin j is 


= I — 2cij^ X uj{ sin^ 




sin2(A|U) 

cos2(f) ^ 


(B48) 


where 


cos(a-^) 




(B49) 


is the effective NV spin rotation angle during one cycle. 

Here the vectors represent the sample nuclear 

spin Hamiltonians in the two subspaces of the NV elec¬ 
tronic spin, i.e., i takes the value of the NV spin state -I, 
0, or I. In the case of nuclear spin-1/2, we have Wg = 
where wl is the nuclear spin Larmor frequency. On the 
other hand, -I- A/, where A/ is the dipolar cou¬ 

pling component along the NV z axis. Then the dip in 
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the signal, = 1 — , can be related to contrast by: 

C{t) = Y[S^ =Y[[1-V^] 


n 


(B50) 


The expression can be further simplified in the limit 
ojl ^ l^il, where A{ = x 4 |[cos(/ 3 sind,sin(psini?,cosd]. 
Then, to second order in A{, the signal is determined by: 

57 ^ sin"(d) sin^ (‘^) sin^ (^) ^ 


COS' 


(^) 


For simplicity in the following steps, we define Kj = 
x4-isin(#) = (Ai^)^ + (Aiy)'^. We can also simplify Eq. 
(B35) using all k values to get 


|5(wL,r)|2 = 


16 sin^ (^) sin" (^) 


Ur 


COS' 


(^) 


The term Kj can converted into an integral of the 
form f p(r^K^(r)d^r and integrated over the sample. 
Since A^ represents the frequency shift from dipolar cou¬ 
pling, one can show from the definition of k that: 

%-fl ^ (1 - (ulf) = ^-feBius- 

3 3 

(B58) 

This along with the approximated expression of the filter 
function finally allows Eq. (B57l to be written as: 


C'(r) « exp 
= exp 




- -"fl{NTf‘SlTLC‘ 

TT^ 




B 


RMS 


(B59) 


Importantly, the expression (B59l for contrast exactly 
matches that given in Eq. (B371, showing the equivalence 


(B52) of the two calculational approaches presented here. 


Then the NV signal contrast from an ensemble of nuclear 
spins precessing at Larmor frequency ujl is 

Cir) = Ylj (1 - il5(wL,r)pK2) . (B53) 

This product can be reconciled with the exponential form 
of the previous section in the following manner. Eirst a 
variance of the effective field is defined as: 

{n^) = -Y^Kl (B54) 

i=i 

The variance is just an average of the individual val¬ 
ues. If the number of nuclear spins n is large, one can 
assume that each spin acts like an average spin, and 
can be replaced with (k^). Then the product simplifies 
to: 


C{t) = rij (1 - l\9i^L,T)\‘^Kj) 




Substitution with Eq. (B54) yields: 


C{t) = I 1 - 


(B55) 


(B56) 


Appendix C: Estimated Proton Nanoscale NMR 
Linewidth Calculated from Correlation Time 


The NV NMR protocol detects a nuclear spin signal via 
the dipole-dipole interaction, which makes it extremely 
sensitive to changes in nuclear spin position. As a conse¬ 
quence of the strong distance dependence of dipolar cou¬ 
pling, nuclei diffusing in a liquid on the diamond surface 
move in and out of the nanoscale sensing volume very 
quickly, which limits the interaction time between the 
NV and nuclear spin. As a result, the nanoscale NMR 
linewidth is broadened. This is in contrast to conven¬ 
tional NMR detection via an inductive coil surrounding 
the sample, in which the nuclei can be fully contained 
within the sensing volume and changes in nuclear posi¬ 
tion have little effect on the signal. 

We assume that the interaction between the NV and 
nuclear spin lasts for a characteristic correlation time, t^, 
and that the probability of finding the particles interact¬ 
ing drops off exponentially in time. By taking the Fourier 
transform, this behavior produces a Lorentzian lineshape 
L(ui) typically written as: 


L(w) 


1 Tg 
TT 1 -k 


(Cl) 


Note that for large n this is the definition of the expo¬ 
nential. Then 


C'(t) = lim I 1 - 

r>—vnn \ g Jl < ^ 


= exp [-i| 5 (a;L,r)p^/ 


(B57) 


This can also be written in a standard Lorentzian form: 


L(iJ) 


1 l/Td 
7ra;2 -h l/rj' 


(C2) 


The full width at half maximum (FWHM) is then 2/rd. 

The translational diffusion correlation time for two 
spins in three dimensions (in our case the immobile NV 
and diffusing nuclei in molecules in the sample) can be 
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FIG. C.l: Estimate of proton nanoscale NMR linewidth as 
a function of NV depth, for immersion oil on the diamond 
surface. 


the sample carrying the nuclear spins. 

Low-fluorescent immersion oil is typically composed 
of liquid polybutadiene mixed with smaller amounts of 
paraffins and carboxylic acid ester^^. In one example of 
an immersion oil with kinematic viscosity v = 450 cSiP^, 
the polybutadiene component has an average molecular 
weight of 1600 g/mol. The hydrodynamic radius of the 
molecule is on the order of r ^ 1 nnP^, and the density 
is p ^ 0.9 g/mL. The dynamic viscosity is then 

r^ = pv = 0.405 cP. (C5) 

We use this viscosity as an approximation for the similar 
immersion oil employed in our experiment. Using the 
Stokes-Einstein relationship 


related to molecular geometries and diffusion coefficients 

bjUMZl 


D = 


keT 

67rpr 


(C6) 


Td 



(C3) 


where d is the distance of closest approach between the 
two spins and D^v is the average of the diffusion coeffi¬ 
cients for the two spins. Since the NV center is immobile, 
we can assume that its diffusion coefficient is zero. The 
distance of closest approach is the NV depth, dNV- Then 
the correlation time becomes: 


Td 


0p2 

k^nuc 


(C4) 


where Dnuc is the diffusion coefficient of the molecules in 


gives a diffusion coefficient Dou « 5 x 10“^^ w?/s. 

Figure C.l plots the estimated nanoscale NMR 
linewidth for immersion oil as a function of NV center 
depth calculated using equation |C4[ The estimated NMR 
linewidth is ~ 5 kHz for a ~ 10 nm deep NV center, while 
the broadest NMR linewidth we expect to see in the mea¬ 
surements performed in this work is ~ 30 kHz for a ~ 4 
nm deep NV center. Consequently, we expect that the 
NV NMR detection bandwidth is much broader than the 
sample’s NMR linewidth (i.e., the infinite approxi¬ 
mation is valid) for nearly every measurement, excepting 
measurements with long pulse sequence durations on the 
shallowest NV centers. 
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